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Abstract: We describe a novel method of tracking the rotational motion 
C^ . of clusters of colloidal particles. Our method utilizes rigid body transfor- 

mations to determine the rotations of a cluster and extends conventional 
proven particle tracking techniques in a simple way, thus facilitating 
cj^ ' the study of rotational dynamics in systems containing or composed of 

O . colloidal clusters. We test our method by measuring dynamical properties 

of simulated Brownian clusters under conditions relevant to microscopy 
^ . experiments. We then use the technique to track and describe the motions 

of a real colloidal cluster imaged with confocal microscopy. 
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1. Introduction 

Suspensions of spherical colloidal particles have proven a valuable model system for under- 
standing many complex phenomena. Perhaps most notably, the model hasprovided insight into 
dynamical processes within different phases of matter j l l I2I j j, 1 4 [sl, lo, LZl, m , as well as various 
mechanisms involved during phase transitions 191 [lOl [ill I12I [ij, ll4J] . A combination of digital 



video microscopy and computerized particle tracking algorithms |tl5l UM has allowed for the 
direct visualization of such colloidal systems and measurement of static and dynamical proper- 
ties under many experimental conditions lll7[ll8[ll9[l2Ql l2in. Given the spherical symmetry of 
the particles, most previous studies have understandably focused only on understanding transla- 
tional dynamics, though there have been a few studies on rotational dynamics in dilute systems 



of spheres 1*221 12311. and on the translational and rotational dynamics of anisotropic colloidal 
particles [24, 2il2il27|]. 



Within the past decade, researchers have developed a variety of techniques for synthesizing 
clusters of colloidal particles with a wide range of reproducible morphologies 1281 129|], several 
of which are shown in Fig.[T] Colloidal clusters have the potential to extend the colloidal model 
past one of simple spherical particles and into a realm where the collective dynamics of particles 
with complex shapes, more representative of molecules, can be studied 11301 1311 13211. 
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Fig. 1. Volumetric images of colloidal clusters with n = 4, 5, and 6 from 3D confocal 
micrographs. Images have been filtered and enhanced to allow easy visualization of the 
3 -dimensional structures. The cases n = 4 and ;i = 5 are accurate representations of the 
simulated tetrahedra and pentahedra discussed in the text. Individual particles are approxi- 
mately 2 /im in diameter. 

Conventional particle tracking methods ||l5l ll6|l are designed to follow the translational mo- 
tions of individual particles, and so are immediately applicable to colloidal clusters, provided 
the particles can be reliably distinguished. However, a description of dynamics within systems 
composed of or including clusters is incomplete without knowledge of how rotational degrees 
of freedom are explored. Such an understanding could provide further insight into fundamental 
behaviors of systems with orientational order, such as liquid crystals, or systems subjected to 
external fields and anisotropic flows ll24ll . To our knowledge, there has been only one study that 
combined video microscopy with particle tracking to measure the rotational motion of colloidal 
clusters 13311, using different methods than described in this work. That study focused on dilute 
systems of planar clusters undergoing two-dimensional diffusion near a boundary, and revealed 
a decoupling of translational and rotational diffusion due to hydrodynamic effects. 

We present here a simple and generalized method to track and calculate the two- or three- 
dimensional rotational motions of clusters of colloidal particles. Our method uses existing par- 
ticle tracking routines and rigid body transformations to measure the changes in orientation of 
a cluster over time. To demonstrate the effectiveness and accuracy of the method, we simu- 
late the motion of colloidal clusters under conditions that are relevant to common microscopy 
experiments. We then use our method to measure the rotational diffusion coefficient of a real 
colloidal cluster, and compare the results of both simulations and experiments to other methods 



of measuring rotational diffusion. 

2. Calculating Rotations 

Our procedure for calculating rotational displacements of clusters is based on a method by 
Challis for determining rigid body transformations between reference frames [34]. Challis' 
procedure was originally intended for comparing osteometric measurements in biomechanical 
analyses. As we will show, the method can be intuitively adapted to measure rotational dynam- 
ics. First, we give a brief reprise of Challis' method, and then describe how it is used to study 
the systems mentioned here. 

2. 1 . Challis ' Procedure for Coordinate Transformations 

Given a set of points which have coordinates {x/} measured in one reference frame and coor- 
dinates {y;} measured in a second frame, there exists a transformation 

y,- = ^Rx,- + v, (1) 

where subscript / refers to the /-th point in the set, ^ is a scale factor, R is a 3 x 3 rotation matrix, 
and V is the vector separation of the two reference frames. For our purposes, we may set the 
scale factor to unity and assume that both coordinate frames share a common origin, thereby 
setting all elements of v to zero. 
We are therefore left with 

y.- = Rx,-, (2) 

and only the rotation matrix to describe the transformation between coordinate frames. The 
rotation matrix is an orthonormal matrix with the properties 

RR^ = RR^ = R^R = E, (3) 



det(R) = +l, (4) 

where E is the identity matrix and det () denotes the determinant. 

For a set of n points, R can be calculated using a least squares approach. This method mini- 
mizes the quantity 

iXl[y,--Rx,f [y,--Rx,-]. (5) 

Ignoring the factor of \/n, expansion of Eq. (O yields 

i,(yifyi + {^if^i - 2(yO^Rx.-. (6) 

Given that all x/ and y^ are fixed, minimizing Eq. Q is therefore equivalent to maximizing 

X:(yO^Rx,- = Tr ( R^- i:yKxO^ ) = Tr(RTc), (7) 

i=\ \ ^ i=\ J 

where C is the cross-dispersion matrix calculated from 

C = j^yi{xi)'. (8) 

i=\ 



At this point, a singular value decomposition is performed on C such that 

C = UWV'^, (9) 

where W is a diagonal matrix containing the singular values of C, and U and V are orthogonal 
matrices. As was shown in ll34ll . upon maximizing Eq. (|7]), R is given by 

"10 1 

R = U 1 V'^. (10) 
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This procedure is applicable to all non-colinear sets of points with n>3. 

2.2. Application to Colloidal Clusters 

Particle tracking [[isl UM yields vector coordinates x • for every particle / over a distinct set of 
times. Hence, the first step in tracking rotational motion of a cluster is to track the translational 
motion of each particle within the cluster. For each cluster, we first determine the center of mass 
xcM at a given time and subtract this quantity from the coordinates of particles belonging to the 
cluster, thereby removing any translational motion. We are left with new coordinates x/ in the 
center of mass frame, 

X,=X--XCM. (11) 

This step is equivalent to setting the elements of v to zero in Eq. ([T]). With translational motion 
removed, we may apply Eq. (O with a slightly different interpretation. Rather than representing 
a transformation between coordinate frames, we may understand R as describing the rotational 
trajectory of a particle with inital position x^ to a final position X; such that 

Xi = Rxl (12) 

Therefore, we may use Challis' procedure to calculate a rotation matrix for each pair of suc- 
cessive times [t^t-\-At]. With the complete set of rotation matrices, {R^}, we may reconstruct 
the entire trajectory of a particle about the cluster center of mass by computing the product of 
successive rotations. Given x^, the position of a particle at some later time t can be calculated 
as 

xi{t) = R^-^^R^-2^^.. rV = n^^x^ (13) 

k 

where the index k enumerates the rotation between successive times. 

The advantage of calculating {R*^} is that it describes the collective behaviors of particles 
within a cluster, rather than a property of any individual particle. For example, knowledge of 
{R"^} for a cluster allows for immediate calculation of the motions of any particular particle 
about the center of mass, or the motion of the cluster about any arbitrary axis of rotation. 
Diffusive anisotropic clusters with large aspect ratios rotate more freely about a long axis than 
about a short axis. Given {R*^}, however, one needs only the initial orientation of these axes to 
compute and compare the motions around them. 

3. Tests of the Prescribed Method 

For simplicity and to better reproduce the data collection process in typical microscopy experi- 
ments, we adopt a length scale of microns and a time scale measured in timesteps (ts), which is 
equivalent in microscopy to video frames or image stacks. 

In conventional particle tracking experiments, it is important to minimize and understand 
the uncertainty, i.e. the noise, inherent in locating a particle. Typical microscopy experiments 



combine high magnification optics with CCD cameras to record raw digital images. In the 
absence of other sources of noise, the uncertainty in particle position, i.e. the minimum noise 
level, depends on the optical resolution [pixels/distance] of the instrumentation and the size 
[pixels] of the object being tracked. Standard image processing and particle tracking techniques 
can locate the centers of particles to within ^ 1/7V of a pixel, where N is the width of the object 
in pixels. Optical resolution varies between experimental set-ups but is typically in the range of 
0.2 /im/pixel. The minimum uncertainty in particle position is the product of these factors. For 
example, observing a 10 pixel wide object with an optical resolution of 0.2 /im/pixel leads to 
a lower limit of ^ 20 nm uncertainty in particle position. Other sources of noise, such as stray 
light entering the microscope, noise within the CCD camera itself, etc., slightly increase the 
uncertainty in particle position and further limit particle tracking resolution. 

To test our method, we simulate the rotational Brownian motion of tetra- and pentahedral 
clusters with different rotational diffusion coefficients, Dr, and different levels of noise, CJ^. We 
first generate noise-free cluster trajectories. For tetrehedra, we place particles at initial coordi- 
nates 

{R/V3)-{l,±l,±l), 

(7?/V3)-(-l,±l,Tl), ^ ' 

where R is the distance from a particle center to the cluster center of mass. For pentahedra, we 
use initial coordinates 

^•(1,0,0), 

(7?/2).(-l,±V3,0), (15) 

7?-(0,0,±l). 

Once initialized, we evolve each simulation for 10^ time steps. 

At each time step, we select three random angles. Of, j3, 7, from a Gaussian distribution with 
a standard deviation of \/2Dr. This distribution ensures that the simulated dynamics will be 
in agreement with the Stokes-Einstein-Debye relation discussed later [Eq. ([2Qb l . Each particle 
in a tetrahedron is rotated by an angle 7 about the z-axis, then by an angle j3 about the j-axis, 
and finally by angle a about the x-axis to produce the tetrahedron at the subsequent time. The 
rotation matrices used are 



R^(^a)=\ cosa -sina , Ry{^) 






Rz{y)= sin7 C0S7 . (16) 



Thus, given an initial position vector x^, the subsequent position vector is x = RxR3;R2X^. 

Trajectories generated in this way are noise-free in the sense that they are absent of un- 
certainty in particle position to within machine precision. To mimic the type of experimental 
noise previously mentioned, we post-process the trajectories by adding Gaussian random num- 
bers, with standard deviation CTx, to the particle coordinates. The levels of noise presented here 
correspond to uncertainties of G^ G {10, 30, 50, 70, 100} nm in the x-, y-, and z-directions. 
Experimental uncertainties are typically within the range of 20-60 nm, and so the levels studied 
here are relevant to microscopy experiments. After the noise is added, we apply our method of 
measuring rotational motion in order to gauge the effect of experimental noise on results. 



4. Analysis of Rotational Motion 

In this section, we calculate the rotational motions of simulated colloidal clusters using the ma- 
trix methods described above. We focus on a method that uses the rotation matrices to determine 
the motion of a fictional orientation vector attached to the cluster. The rotational displacements 
of one or more such vectors about the center of mass can be used to compare and contrast mo- 
tions about different axes, which is especially useful when the objects studied are anisotropic 



and motions about axes are expected to differ (SJ, |3 



We note that a second method exists to calculate rotational displacements using solely the 
rotation matrix. This method determines the rotational axis from R and can then be used to 
calculate the magnitude of an angular displacement. While a mathematically direct and general 
approach, it describes only the average cluster dynamics and can provide no insight into how 
motions about different axes vary. Given these limitations, we relegate a description of this 
method to the appendix. 

4.1. Rotations via Orientations 

To analyze the rotations of a cluster, we consider the motions of an initial orientational unit 
vector p^ fixed to the cluster. We determine its orientation at a later time t by applying the set 
of rotations such that 

p(o=nRV, (17) 

k 

in a similar way as in Eq. ([T3l) . Although there are no constraints on what one may select p^ to 
be, some choices may be more enlightening than others. For example, the cluster n = 5 shown 
in Fig. [T] has distinct long and short axes, and so one expects slower diffusion about the short 
axis. To quantify how rotational dynamics about these axes differ, one could choose two p^ to 
study separately: one choice of p^ perpendicular to the long axis; a second perpendicular to the 
short. Such a procedure would yield information relating to motions about the long and short 
axes, respectively. 

In line with Il35l|36l|37|], we define a vector rotational displacement 

cp(0 = / a){t')dt' (18) 

Jo 

in the time interval [0,^]. The vector (b{t')dt' has a direction given by p(^') x ^(t' ^ dt') and 
magnitude \a){t')dt'\ = cos~^ [p(^0 • ^{t' ^dt')], which is the angle subtended by p during this 
time interval. 

To illustrate the meaning of $(0. consider an object with constant angular velocity O) = (oz. 
Over a time A^ the rotational vector displacement is given by $(Ar) = a)^t = (0,0, coA^). 
Therefore in general, (p(^) has components in each of the Cartesian axes, ((p;c, ^>y^ ^z), corre- 
sponding to cumulative rotations about those axes. 

In Fig.Oa), we show the orientational trajectories of particles within a simulated tetrahedron 
projected onto the surface of a unit sphere. Fig. Ob) shows 2D projections of trajectories for 
two p^ through the rotation space described above. As shown in Fig.[2l even though the cluster 
is a solid body, trajectories of individual particles differ due to rotations of the cluster about 
random axes. 

Given the definition of (p(^) in Eq. ([TSl) , we may define an unbounded mean square angular 
displacement (MS AD), akin to a translational mean square displacement, as 

(A(p2(A0) = ([cp(^ + A0-$(0]'), (19) 




Fig. 2. (a) Trajectories of particles within a simulated Brownian tetrahedral cluster pro- 
jected onto the surface of a unit sphere. Colors represent different particles, i.e. different 
choices for p^. (b) 2d projection of two trajectories through rotation space {(px,<Py)- Colors 
correspond to the same particles in (a). Both trajectories in (b) begin at (0,0) and end at 
open circles. 



where the angle brackets indicate an average over all equivalent lag times A^ In three dimen- 
sions the Stokes-Einstein-Debye relation states that the MS AD grows as 



{A(p^{At))=4DRAt, 



(20) 



where Dr is the rotational diffusion coefficient. In this paper, we focus on measuring the 
MSADs of clusters, but we point out that other techniques exist to quantify rotational dynamics. 
For example, by observing the decay of an orientational correlation function (p(^ + A^) • p(^)) 
one can measure Dr 1381 |39|]. Our method can be applied, in this case, to compute p(^) as in 
Eq. dnl). 

In Fig. Oa) & (b), we show the MS AD of two simulated tetrahedral clusters with diffusion 
constants 10~^ rad^/ts and 10~^ rad^/ts respectively and different levels of noise. The influence 
of noise is apparent in Fig. [3] as deviations from linearity at small A^ Eventually, the MS AD 
recovers the true diffusive behavior because the cluster has made rotations large enough to 
distinguish from the noise. 

With no experimental noise, a log-log plot of MSAD as a function of At will be a straight 
line with a slope of unity, as indicated by the open circles in Figs.Oa) & (b). However, in the 
presence of noise in particle positions, rotations cannot be accurately resolved below a certain 
threshold, <I>. For example, a stationary cluster will appear to make small, but fictional, rotations 
as a result of this noise, and measurements of the MSAD will yield 



(Acp^(AO)=<I>^ 



(21) 



where <I>^ is independent of At due of the lack of correlations in noise. 

In the case of translational diffusion, one expects that the noise in each direction will con- 
tribute an error of 2a^ to the mean- square displacement llldl . In terms of an angular uncertainty, 
this contribution is diminished by a factor ofR^, where R is the average distance from a particle 
to the cluster center of mass. Our matrix method reduces this uncertainty further by a factor of 
n. Thus, when only static noise is present, one expects for a diffusing cluster 



{A(p^{At))=4DRAt^^^ 



(22) 
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Fig. 3. Mean square angular displacements of simulated tetrahedral clusters for different 
noise levels Ox and diffusion coefficients of (a) 10~^ rad^/ts and (b) 10~^ rad^/ts. Open 
circles are the theoretical MS AD based on Eq. ([20l) . Deviations from linearity at small 
A^ demonstrate the effect of noise when resolving small rotations. Deviations at large Af , 
however, are the result of low statistics at these lag times. 



where 



^^ = ea^/nF?- 



(23) 



To test this assertion, we add noise to simulations of stationary (Dr = 0) tetrahedral clusters 
and calculate the MSADs. As in Eq. (|2T1) . the MSADs are constant over time. We take <l>^ to be 
the value of the MS AD as At -^ 0. Shown in Eig.jlja), Eq. (1231) accurately describes the static 
angular uncertainty for a wide range of tetrahedral cluster sizes, R, and noise levels. Fig.Ub) 
shows all of the previously measured MSADs in Fig. [3] plotted with the noise subtracted. This 
precisely collapses the MSADs to the true values in each case. Values of <J>^ for pentahedral 
clusters are also predicted by Eq. ([23l) , and the MSADs for these clusters can be similarly 
collapsed. 

Knowing how uncertainty in particle positions affects measurements of dynamics quantities, 
such as the MS AD, is clearly important. However, one is typically unable to define the noise 
level so precisely in experiments. As stated earlier, particle tracking resolution depends on var- 
ious experimental factors including particle size, optics, and the type of camera used to record 
images. As shown in Fig.Oa), for slowly diffusing clusters and moderate-to-high noise levels, 
the measurements of the MS AD may not be representative of the true dynamics until fairly large 
lag times. At these large At, statistics are poorer than small At, and confident measurements are 
harder to obtain. Therefore, this type of oversampling can be detrimental to microscopy exper- 
iments where data sets consist of a few hundred to a few thousand images. Thus it is important 
to determine an imaging rate that yields the most information from a data set. 
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Fig. 4. (a) Measured values of 4>^ for non-diffusing tetrahedral clusters. Colors indicate 
noise levels and symbols indicate cluster size. Solid line is the prediction <I>^ = 6o^ /nR^. 
(b) Same data in Fig. [3] with the appropriate O^ for each noise level subtracted. ^^ for 
pentahedral clusters collapses in the same way. 



In microscopy experiments, one can effectively increase the rotational signal-to-noise ratio 
between each time step by imaging less frequently. Determining the appropriate sampling rates 
in experiments can be done by estimating the noise level cr^, the diffusion constant Dr and cal- 
culating an approximate <I>^. Diffusive motion will begin exceeding noise when <l>^ '^ ADr^I. 
Solving for A^ then yields a reasonable sampling lag time. It is also important to avoid under- 
sampling, that is, too long a lag time. Doing so will make diffusive motion between images 
appear erroneously slow. We find that a good rule of thumb for the upper limit on sampling 
time should be the time when the cluster has diffused ^1 radian^. Thus, an estimate for the 
upper limit is ^ \/ADr. However, as in all particle tracking experiments, the time between im- 
ages must be small enough that individual particles can be confidently identified. This typically 
means that particles must be imaged before moving a distance of one interparticle spacing [16]. 

5. Experimental Application 

Measurements of rotational motion of diffusive tetrahedral clusters have been performed using 
the described method 112 Sll . Real fluorescent tetrahedral clusters are synthesized as in 112 Sll . A 
cluster is composed of individual poly(methyl methacrylate) (PMMA) spheres, each with a 
diameter 2.45 /im as measured by static light scattering (SLS). The particles within a cluster are 
irreversibly bound together, but are sterically stabilized to prevent the possibility of aggregation 
to other clusters. Dilute suspensions of clusters are prepared in a mixture of cyclohexyl bromide 
(CXB) and c/^'-decalin (DCL) at a ratio of 85/15 (w/w) that closely matches both the density 
and index of refraction of the particles. Clusters are imaged in 3D over time with a Leica TCS 
SP5 confocal microscope. We track locations of the individual particles within a tetrahedron 
using standard particle tracking routines Ill5lll6|]. The uncertainties in particle position for these 
experiments are ?^ 30 nm in the x- and }^-directions, and ^ 40 nm in the z-direction. Given these 
tracking resolutions and assuming a maximally packed tetrahedron, from Eq. (l23l) we estimate 
the angular resolution in this experiment as ^ ?^ 0.028 radians (1.6°). 

Once tracked, we calculate the translational MSD and the MSAD and determine the trans- 
lational and rotational diffusion coefficients, Dj and Dr, respectively. In three dimensions, the 
translational MSD is described by the Stokes-Einstein- Sutherland equation 
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Fig. 5. Calculated MSD (blue squares) and MS AD (red circles) of a diffusing tetrahedral 
cluster. Solid lines are best fit lines over range of data used to determine the respective 
diffusion coefficients. 



(Ar^) = eDrAt, (24) 

while the MS AD is described by Eq. (l2Ql) . 

Hydrodynamically, tetrahedral clusters can be accurately modeled as spheres 140"] given the 
relation 



<^tetra = 1-844 X ^sphere, 



(25) 



where Jtetra is the effective hydrodynamic diameter of the cluster, and ^sphere is the diameter of 
the particles within the cluster. Theoretical translational and rotational diffusion coefficients, 
Dt and Dr respectively, can be calculated using 



Dt = 



Dr 



kBT 



3717] dtt 

ksT 



(26) 



(27) 



^^4tra' 

where Rb is Boltzmann's constant, T is the temperature, and T] is the viscosity of the suspending 
solvent. The viscosity of the CXB/DCL mixture was measured at T] =2.18 mPas and experi- 
ments were performed at T = 295 K. 

In Fig.[5l we show the MSD and MS AD of a diffusing tetradedral cluster. The MSD corre- 
sponds to motions of the center of mass of the cluster, and the MSAD is an average over the 
individual particle MSADs (i.e., using the initial orientations of the particles as four separate 
p^). Solid lines are fits to the data over the indicated range 3 s < At < 22 s . As can be seen, 
both the MSD and MSAD are approximately linear only at small lag times. This is because the 
data set used to make these measurements consists of only a single cluster for less than 300 
timesteps (^ 220 s). Diffusion coefficients are measured by performing linear fits to the MSD 
and MSAD shown in Fig. [5] over the indicated range and extracting the slopes of the lines. 

Shown in Table [T] are the measured translational and rotational diffusion coefficients. Values 
of dtQtra are calculated from Eqs. (l26l) & (l27l) using known experimental conditions, and values 



Table 1 . Measurements from Tracking Tetrahedral Cluster 
Diffusion Coefficient Jtetra (M^^) <^sphere (M^^) 

Trans. (3.85 ±0.49) x 10" V^^/s 5.15±0.67 2.79±0.13 

Rot. (4.93 ± 0.49) x 10"^ rad^/s 4.94 ± 0.20 2.68 ±0.11 



for ^sphere then follow trivially using Eq. (1251) . Ideally, the sizes calculated from translational 
motions would be identical to those calculated from rotational motions. As shown, these values 
agree to within 4%. Such good agreement between these two measurements demonstrates our 
ability to track translational and rotational motions of clusters simultaneously. 

Previous work using standard particle tracking and our matrix method reported particle diam- 
eters ^ 14% larger than the 2.45 /im measured by SLS 112 811 . A reanalysis of the data, presented 
above, improves the particle diameter measured using our matrix method to within ^9%. We 
attribute the better accuracy of using rotational motion, as opposed to translational, to the least 
squares minimization process, which incorporates the rotational motions of all particles within 
a cluster. The remaining disagreement is likely due to multiple factors, including the fitting rou- 
tines used in SLS and swelling of particles when in a solvent of CXB/DCL. SLS measurements 
were performed with suspensions in pure DCL, which is not known to swell particles. 

6. Conclusion 

We have presented a simple method of tracking the rotational motions of colloidal clusters. 
Our method implements conventional particle tracking routines to determine the locations of 
individual particles within a cluster and uses this information to compute rigid body transfor- 
mations that describe changes in a cluster's orientation over time. The set of matrix transforma- 
tions constitutes a global description of a cluster's motion during the course of an experiment or 
simulation and allows one to calculate rotational dynamics about any arbitrary axis. The least 
squares minimization used in this method considers the motions of each individual particle 
in calculating the rotation matrix, and therefore measurements of rotational displacements are 
less sensitive to tracking noise. Additionally, precision in determining angular displacements 
increases with cluster size R. When combined, the resolution of measuring angular displace- 
ments scales with llnR^. 

Though constructed specifically for tracking colloidal clusters, we emphasize that our 
method is not limited to these systems, but is applicable to following the rotational motion 
of any body over time, providing that at least three distinct noncolinear points in the body can 
be reliably distinguished. Because the accuracy of tracking a rotation depends on the number of 
particles and the size of the body, the possibility exists that this method can be adapted to many- 
body systems where each particle diffuses independently of the rest while the system itself also 
undergoes bulk rotations. In this scenario, diffusive motions are random and can therefore be 
treated as noise while calculating {R^} for the system. The rotation matrices can be inverted 
and, in a manner similar to Eq. ([T3l) . the bulk rotations can be removed, leaving only uncorre- 
lated diffusive motion. Thus, this method is also applicable to determining particle motions in 
rotating coordinate frames. 
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Appendix 

Given a rotation matrix R, one can calculate at each time the axis of rotation u and angular 
displacement Acp. In this notation, a rotational displacement can be described by a vector A(p = 
A(pu, where u has components in each of the Cartesian axes. 

A single rotation will, by definition, have no effect on the direction of u, therefore. 



Ru = u. 



(28) 



From Eq. (1281) . we see that the axis of rotation is an eigenvector of the matrix R with an eigen- 
value of 1. For a set of rotational displacements {R^}, one may determine the axes of rotation 
by calculating the eigenvectors and eigenvalues of the rotation matrices, searching for the eigen- 
values equal to 1, and taking the corresponding eigenvectors. 

To determine the size of the displacement about the axis of rotation, one defines an arbitrary 
vector w perpendicular to u. For simplicity, we choose w to be perpendicular to the x-axis 
(denoted by i). 



W: 



UX ] 



UX ] 



and apply the rotation matrix, 



Rw = w'. 



(29) 



(30) 



The magnitude of the displacement is the angle between w and w', and can be computed using 
the cross product relation 



sin (A(p)u = w X w'. 



(31) 



Calculating displacements relative to the axis of rotation always results in displacements 
greater than or equal to those measured relative to an arbitrary p^. For example, measuring a 
diffusion coefficient of a spherically symmetric body with this method will yield a value that is 
a factor of 3/2 of the actual diffusion coefficient in Eq. i 
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Fig. 6. The MSADs of pentahedra with diffusion coefficients Dr = 10~^ rad^/ts and Dr = 
10~^ rad^/ts and Gx = 50 nm. Dashed lines show the uncorrected MS AD described below 
and solid lines are corrected data, obtained by multiplying the dashed lines by 2/3. Open 
circles are the theoretical MS AD from Eq. ( |20t . 



To understand the origin of this difference, consider the rotational axis u = z, a perpendicular 
vector w = X, and an arbitrary vector p^ located at spherical coordinates (sin 0, 0, cos 0). If w is 



rotated by an amount Scp^, the angle between w and w' is identically 5(pu. However, the angle 
5 (p between p and p' can be shown to be 

5cp = cos~^ (sin^0cos5(pt^ + cos^0) (32) 

For small 5(pt^, we can approximate 5(p as 

1 - 5(p^ ^ sin^ e (1 - 5(pI) + (1 - sin^ O) (33) 

5(p^^sm^e5(pl. (34) 

If an average of Eq. (l34l) is taken over spherical coordinates, we are left with an expression 
similar to an MSAD, 

{5<p') = {sm^d){5<pl) = ^{5<pl), (35) 

The factor of 2/3 arises for the same reasons in Perrin's original derivation of rotational 
diffusion 114 IH . Fig. [6l shows the MSADs measured in this way for two pentahedral clusters. 
As stated, one can correct for the overestimation of motion by simply multiplying the MSAD 
by 2/3. The average MSAD for any body can be corrected in the same manner, however we 
stress that this method returns only the average dynamics. While more direct than using an 
orientation vector, characterizing anisotropic bodies in this way will convolute motions about 
separate axes, and so will require some care when interpreting results. 



